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SUMMARY 

The use of digital computers in performing non-numeric operations is 
rendered possible by the evolution of new computer languages such as Formac. 
In this paper the problem of using a digital computer to derive the equations 
of motion of a particle in a general orthogonal curvilinear coordinate system 
is considered. Since this operation involves a formulation in terms of first 
and second order differential coefficients, it provides a good demonstration 
of a computer's capability to do non-numeric work, and to assist in the 
formulation process which normally precedes the numerical data processing 
stage. Moreover, this particular problem selves to illustrate the advantages 
of the mathematical techniques employed. Because of the invariant nature of 
the formulation with respect to coordinate transformations, these techniques 
can be used to reduce complicated formulation problems to routine computer 
operations. In applying this procedure, the user need only know the coordi- 
nate transformation equations relating the curvilinear coordinates to an 
orthogonal Cartesian set. A computer program has been written to implement 
these ideas. When this program is ’used and the coordinate transformation 
equations are supplied as input, the computer will output the equations of 
motion. The equations of motion obtained will be relative to the curvilinear 
coordinate system specified by the coordinate trans format ion equations used 
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as input. Results are presented for the following curvilinear coordinate 
systems: spherical polar, cylindrical polar, oblate spheroidal and prolate 
spheroidal. 

HffRODU CTI ON 

Research undertaken with the object of promoting man-computer inter- 
action has directed attention to the use of computers for non-numeric 
operations. In particular, the possibility of using digital computers to 
derive the equations of motion and of mathematical physics in a general 
curvilinear coordinate system has been explored. Traditionally, these 
functions were considered to be the exclusive preserve of the scientist, 
nevertheless, as is shown in this paper, digital computers can participate 
in the performance of such tasks. However, if the extensive logic and stor- 
age capabilities of these computers are to be used to full advantage, a 
departure from conventional techniques of formulation may be necessary. The 
extent to which conventional methods should be modified to enable digital 
computers to participate effectively in non-numeric operations has been 
examined. For example, when conventional methods are used, the form which the 
equations of motion and of mathematical physics assumes depends on the 
coordinate system used to describe the problem. This dependence, which is 
due to the practice of expressing vectors in terms of their physical compo- 
nents, can be removed by the simple expedient of expressing all vectors in 
terms of their tensor components. As a consequence of the geometrical 
simplification inherent in the tensor method, the operations involved in 
formulating problems in unfamiliar curvilinear coordinate systems can be 
reduced to routine computer operations. It is this aspect of the tensor 
method which makes it so attractive for the types of computer applications 
contemplated in this paper. 
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NOMSNCLATURE 
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d 2 xi 


dt 
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time 
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system coordinates 
physical component of force 
potential function 
gradient of potential function 

Super cripts 

indices of contravariance 


ANALYSIS 

A formulation of the equations of motion of a point mass, vhich is valid 
in all orthogonal curvilinear coordinate systems, may be obtained by expressing 
all relevant vectors in terms of their tensor components, rather than in terms 
of their physical components. This formulation gives rise to the tensor 
equations of motion. In order to indicate the method of dealing with mixed 
quantities, part of the force system is assumed to be given in the form of 
the gradient of a potential function. Such a force assumes the covariant 
form (ref. l) . The existence of forces which are known only in terms of 
their physical components is also assumed. Moreover, since acceleration and 
velocity are contravariant vectors, the equations of motion have to be 
formulated from a system of covariant, contravariant and physical quantities. 
Since compatibility requires that the two sides of every equation must bal- 
ance with respect to their covariant or contravariant properties, it is 
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necessary to convert all the force terms to the eontravariant form. When 
appropriate conversions are made, the resulting equations are in a form well 
suited to non-numerie computer operations. 

Omitting the details of the derivation, the procedure may be described 
as follows: let the coordinate transformation equations relating a curvilinear 

coordinate system x, to an orthogonal Cartesian system y, be given by 


y 1 = jrKx^xS) i = 1,2,3 (l) 

In terms of the first and second order partial differential coefficients 

of y with respect to x, the equations of motion of a particle of mass M 
relative to the curvilinear coordinate system x assume the following form: 

JT ’ 
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^ a -n a 

ay oy 


J 


( 2 ) 


where i,j,k,a = 1,2,3 and the summation convention for indices is assumed 
to be operative. That is, if a given index occurs twice in any expression, 
the expression must be summed with respect to that index. An exception to 
this rule occurs when repeated indices are enclosed in parentheses. 

Parentheses around an index imply that the s ’animation convention is to be 
suspended for that index. This means that for each value of "i," equation (2) 
must be summed on a,j, and k. For example, when equation (2) is summed on 
a, it appears as follows: 
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( 3 ) 

The left side of this equation must also he summed on j and k. When each 
ot these indices is permitted to take the values 1^2, 3j in turn^ equation (3) 
assumes the following form: 
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This equation is in a form well suited to routine non-numeric computer opera- 
tions. The large number of terms appearing in equation (4) is due to the 
generality of this equation, which is applicable to any space of three 
dimensions. Moreover, since this equation is applicable to any space of three 
dimensions, it may be permanently stored in the computer. Hence, in order to 
obtain the equations of motion in any system of coordinates, the only informa- 
tion required is the special form of equation (l) relating that system of 
coordinates to the orthogonal Cartesian coordinates y 1 . For example, 
consider a transformation of coordinates specifying the relation between the 
spherical polar coordinates x 1 and the orthogonal Cartesian coordinates y 1 . 
In this case, equation (l) becomes: See sketch (a). 

y 1 =x ! sin x 2 cos x 3 
y 2 =x> sin x 2 sin x 3 
y 3 =x ! cos x 2 



sketch (a) 
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These coordinate transformation equations were supplied, as input to an IBM 
709^ computer, which was programmed for non-numeric operations. When the 
computer was instructed to perform the operations involved in equation (4), 
the following output was obtained in Fortran language. 

COMPUTER OUTPUT 

FOR 1=1 

The expression input for Y(l) is given below. 

X( l)*FMCSIN( X( 2) ) *FMCCOS ( X( 3 ) $ 

FOR I - 2 

The expression input for Y(l) is given below. 

X( 1 ) *FMCSIN( X ( 2 ) ) *FMCSIE( X ( 3 ) ) $ 

FOR I = 3 

The expression input for Y(l) is given below. 

X( l) *FMCCOS(X( 2) )$ 

Equations of Motion 

FOR I = 1 

The eq'oation for I = 1 is given below. 

M*(P(1) -R(2)**2.0*X(l) -R(3)**2.0*X(1)*FMCSIR(X(2))**2.0)$ 

=DPHl( l)+TAU(l)$ 

FOR 1=2 

The equation for I = 2 is given below. 

M*( P( 2 ) *X( 1) **2 . OfR( 1) *R( 2) *X( 1) *2 . 0-R( 3 ) **2 . 0*X( l) **2 . 0*FMCSIE( X( 2) ) *FMCCOS( X( 2) ) ) $ 
=DPHI ( 2 )+TAU( 2 ) *X( 1) $ 

FOR I = 3 

The equation for I = 3 is given below. 

M*(P(3)*X(l)->Ht2.0*FMCSIH(x(2))**2.0+R(l)*R(3)*X(l)*FMCSin(X(2))**2.0*2.CM-R(2)*R(3)* 
X ( 1 ) **2 . 0*FMCSIE( X( 2 ) ) *FMCCOS ( X( 2 ) ) *2 . 0 ) $ 

=DPHI( 3 )+TAU( 3 ) *X( 1) *FMCSIR( X( 2) ) $ 
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JOB ACCOUNTING 


comp/load 


TIME 


TIME 


EXECUTIVE 

TIME 


950.09 


In interpreting these Fortran statements, it must he remembered that: 


R(i) = 


P(i) = 


In terms of conventional mathematical symbolism, these equations assume 
the following form: 


M fed - x 1 (sad) - x 1 (sin x 2 - ^L. + t* 

L at 2 \ d v \ dty J 


M (xi) 


2 d 2 x 2 0 , dx 1 dx 2 / n . 2 „ 2 / dx 3 ^ 2 

— - + 2X 1 (xi) sin x 2 cos x 2 — r— 

at 2 dt at V dty 


+ x 1 ^ 


M (x 1 sin x 2 ) + 2X 1 sin 2 x 2 — — + 2( x i) 


dt dt 


2 sin x 2 cos x £ 


dx 2 dx 3 
dt dt 


+ x 1 sin x 2 t 3 


Because of its generality, equation (4) is applicable in all coordinate 
systems. Therefore, in order to obtain the equations of motion in any other 
coordinate system, all that is required is to supply the computer with the 
appropriate coordinate transformation equations. As a further illustration of 
the procedure involved, consider the equations of motion in a cylindrical 
polar system of coordinates. In this case, the coordinate transformation 


equations are: see sketch (b) . 
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y 1 = x 1 cos x 2 
y 2 = x 1 sin x 2 

y 3 - ^ 


When these coordinate transformation equations were used to evaluate the 
terms of equation (4), the following computer output was obtained. 



sketch (b) 
COMPUTER OUTPUT 


FOR I = 1 

The expression input for Y(l) is given below. 
X(1)*FMCC0S(X(2)) $ 
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FOR 1=2 

The expression input for Y(l) is given below. 
X(l)*FMCSIN(X(2)) $ 

FOR I =3 

The expression input for Y(l) is given below. 

x(3) $ 

EQUATIONS OF MOTION 

The equation for I = 1 is given below. 
M*(P(l)-R(2)**2.0*X(l))$ 

=DPHI(1)+TAU(1)$ 

The equation for I = 2 is given below. 

M*( P( 2 ) *X( 1 ) **2 . Of R( 1) *R( 2 ) *X( l) *2 . 0) $ 

-DPHI ( 2 ) +TAU( 2 ) *X( 1 ) $ 

The equation for I = 3 is given below. 

M*P(3)$ 

=dphi(3)+tau(3)$ 

JOB ACCOUNTING 

comp/load EXECUTIVE 

TIME TIME TIME 

ON MIN. MIN. 

037.78 


1.15 


.14 
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Translating these equations from Fortran language to conventional mathematical 


symbolism yields the following: 


M 


d'^’x 1 x dx 2 dx g 


dt c 


dt dt 


M 




d 2 x g + ^ dxi dxf. 
dt 2 dt dt_ 
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d 2 x 3 " 


Sep 

dx 1 

dp 

dx 2 

dp 


+ T 2 


+ X^-T 2 




( 5 ) 


L dt 2 J dx 3 


+ T" 


Prolate Spheroidal Coordinates 

Another interesting system of orthogonal curvilinear coordinates are the 
prolate spheroidal coordinates. Coordinate surfaces are obtained by rotating 
a family of confocal ellipses and hyperbolae about their major axes. Rotation 
of these conic sections gives rise to a system of prolate spheroids and hyper- 
boloids of two sheets. A family of planes through the axis of rotation 
completes the system of orthogonal surfaces. The curvilinear coordinate 
systems generated by the curves of intersection of these surfaces are useful 
in certain quantum mechanical problems. (Ref. 2 ). The transformation equa- 
tions relating this system of coordinates to the orthogonal Cartesian system 
are as follows : 


y 1 = a sinh x 1 sin x 2 cos x 3 
y 2 - a. sinh x 1 sin x 2 sin x 3 
y 3 = a cosh x 1 cos x 2 

In order to obtain the equations of motion relative to a prolate spheroidal 
system of coordinates , these transformation equations were substituted for 
equation (l) in the computer program. Execute time was 1.63 minutes. 
Omitting the print -out in Fortran language, the equations of motion were 


obtained as follows: 
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M 


a 2 (sin 2 x 2 + sinh 2 x 1 ) — + 2a 2 sin x 2 cos x 2 ^2- 

v - o 


dt" 


dt dt 


+ a 2 sinh x 1 cosh x 1 ^77- ^77- - a 2 sinh x 1 cosh x 1 ^77- ^777- 

dt dt dt dt 


a 2 sin 2 x 2 sinh x 1 cosh x 1 

dt dt 


= a J( sin 2 x 2 + sinh 2 x 1 )t 1 + _J^E. 

dx 1 
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a 2 ( sin 2 x 2 + sinh 2 x 1 ) ^ 3*1 - a 2 sin x 2 cos x 2 dx " L dx ' 

dt 2 


* dxi 

dt dt 


+ 2a 2 sinh x 1 cosh x 1 §Z-L + a 2 

dt dt 


sin x 2 cos x 2 — 

dt dt 


a 2 sin x 2 cos x 2 sinh 2 x 


2 c-inV 2 dx 3 

dt dt 
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sin 2 x 2 + sinh 2 x 1 ) t 2 + 
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J*p 

^x 2 
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a 2 sin 2 x 2 sinh 2 x 1 ■ + 2a 2 si 


dt" 


sin 2 x 2 sinh x 1 cosh x 1 ^777- 

dt dt 


+ 2a 2 sin x 2 cos x 2 sinh 2 xi ^L- 

dt dtj 


= a sm 


in x 2 sinh x 1 t 3 


Sep 

S ^ 3 


Ohlate Spheroidal Coordinates 

If a family of confocal ellipses and hyperbolae are rotated about their 
minor axes, a system of surfaces is generated. These surfaces are the oblate 
spheroids and hyperboloids of one sheet (ref. 3 ). These surfaces, together 
with a family of planes through the axis of rotation, constitute a family of 
orthogonal surfaces. The curvilinear coordinate systems generated by the 
curves of intersection of these surfaces are called oblate spheroidal 
coordinates. Oblate spheroids are sometimes referred to as planetary 
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ellipsoids, because the Earth and the planet Jupiter are approximately of this 
form. The transformation equations relating this system of coordinates to the 
orthogonal Cartesian system are as follows: 

y 1 = a cosh x 1 sin x 2 cos X 3 
y 2 = a cosh x 1 sin x 2 sin x 3 
y 3 = a sinh x 1 cos x 2 


These transformation equations take the place of equation (l) in the computer. 
In this case, the time required to execute the program was I.63 minutes. 
Omitting the print-out in Fortan language, the equations of motion relative 
to a system of oblate spheroidal coordinates were obtained in the following 
form: 
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a 2 (sinh 2 x 1 + cos 2 x 2 ) 


2 ~ 2 ' d2xl + a 2 (sinh xi cosh x 1 ) 


dt" 


dt dt 


- 2a 2 cos x 2 sin x 2 - a 2 sinh x 1 cosh x 1 

dt dt dt dt 


2 , . . . . . 2 2 dxf dx 3 

- a cosh x 1 sinh x 1 sin x — — — 

dt dt 


= a 


^ Jslrih 2 x 1 


+ cos 2 x 2 ) T 1 + 


Sep 

dx 1 
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a 2 ( sinh 2 x 1 + cos 2 x 2 ) + a 2 sin x 2 cos x 2 dx^ dxi 


dt" 

dx 1 dx 2 _ 2 . 


dt dt 


.2 dx 2 dx 2 


+ 2a 2 sinh x 1 cosh x 1 — a 2 sin x 2 cos x' 

dt dt dt dt 


- a 2 cosh 2 x 1 sin x 2 cos x 2 


dx 3 dx 3 
dt dt 


= a J(sirih 2 x 1 + cos 2 x 2 ) t 2 + 


dx 2 
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M 


a" cosh 2 x 1 siri x' 


2aax,„ 2 .,i -ui.22 ax 1 dx 3 

+ 2 a sinh x A cosh x 1 sm x 


a d 2 x 3 


dt= 


dt dt 


+ 2 a 2 cosh 2 x 1 sin x 2 cos x 2 

dt dt 


= a cosh x 1 sin X 2 t 3 + 

dx 3 


CONCLUSIONS 


Digital computers can he 'used to perform a wide range of non-numeric 


operations if they are properly programmed. Research indicates that these 
computers can he used more effectively for this purpose, if all vector 
quantities are expressed in terms of their tensor components rather than in 
terms of their physical components. It is a consequence of the geometrical 
simplification inherent in the tensor method, that the operations involved 
in formulating problems in unfamiliar curvilinear coordinate systems can he 
reduced to routine computer operations. To implement the proposed method, a 
digital computer program has been written to perform a variety of non-numeric 
operations. In order to illustrate the ideas embodied in this report, the 
program has been used to derive the equations of motion of a point mass in 
any coordinate system requested by the user. The results are presented in 
Fortran language. However, for the convenience of readers, the Fortran 
statements are translated to conventional mathematical symbolism. The 
exploitation and extension of these techniques should lead to a substantial 
reduction in the man hours required to formulate and process engineering and 


scientific problems. 


-16- 


REFERENCES 

1. Howard, James C.: Application of Intrinsic Differentiation to Orbital 

Problems Involving Curvilinear Coordinate Systems. NASA TK D-2861, p. 4. 

2 . Margenaw, H. , and Murphy, G. M. : The Mathematics of Physics and Chemistry. 

D. Van Hostrand Co., Inc., pp. 5-7* 

3. Korn, G. A. and Korn, T. M. : Mathematical Handbook for Scientists and 

Engineers. McGraw-Hill Book Co., Inc., 1961, p. 177. 



